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1 Introduction 



Calculation of weak matrix elements of hadrons is one of the central subjects in numerical 
simulation of lattice QCD. An important ingredient in such studies is the value of renormal- 
ization factors that relate operators on the lattice to those in the continuum. In this article 
we present a one-loop evaluation of renormalization factors for bilinear and four-quark oper- 
ators for the Kogut-Susskind(KS) fermion action. The results have been used in our recent 
work for the pion decay constant |[l|] and the — mixing matrixQ. 

Perturbative calculation of renormalization factors for KS fermions has been developed 
in several previous studies |P, ||, |], |^. In particular Daniel and Sheard[§, |^ calculated the 
renormalization factors for KS bilinear operators and a subset of four-quark operators]^ 
in the Feynman gauge. For applications in numerical simulations, however, their results 
need to be extended in several directions. (1) Weak operators for KS fermions are generally 
extended in space and time. The calculation of Daniel and Sheard has been carried out for the 
operators which are made gauge invariant through insertion of gauge link variables between 
quark and anti-quark fields. In recent numerical simulations, however, an alternative method 
of evaluating matrix elements of gauge non-invariant operators without link insertions on 
gauge fixed configurations has been employed0]. In fact, whether the two types of operators 
yield consistent results is the question we have recently addressed ||. Analyzing this 
problem requires the renormalization factors for gauge non-invariant operators as well as 
for gauge invariant ones. We have therefore carried out the calculation for both types of 
operators. (2) Four-quark operators of form O2 = {(11(12) (QsiD with a,b the color indices 
generally mix with those of form Oi = {(11(12) {(13(1 a) under renormalization. We evaluate the 
full renormalization factor for the two sets of operators, while the previous work of Sheard]^ 
listed explicit results only for O2 mixing with itself and with Oi. (3) Renormalization factors 
for lattice operators generally take larger values than those for continuum operators due to 
contribution of gluon tadpoles. Lepage and Mackenzie have argued that the tadpole 
contributions can be removed by a rescaling of quark and gluon fields. We work out the 
renormalization factors for the rescaled operators and examine to what extent their values 
are reduced by rescaling. (4) In addition to the extensions above we have carried out the 
calculation in the general covariant gauge which allows us to check the gauge parameter 
independence of the results for gauge invariant operators. 

For bilinear quark operators a calculation similar to ours has been reported recently 
by Patel and Sharpe|TO|. Our results are in agreement with theirs and also with those of 
Ref. for gauge invariant operators. Patel and Sharpe have extended their calculation to 



four-quark operator s|]rT|. For the gauge non- invariant operators which are relevant for the 
K meson B parameter their values fully agree with our results. We also find agreement 
with the results of Sheard]^ for gauge invariant four-quark operators when a comparison is 
possible. 
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We should mention that we do not treat penguin operators in this article. Calculation 
of their renormalization factors is technically feasible, which should be pursued in future 
investigations. 

This paper is organized as follows. In Sec. 2 KS operators whose renormalization factors 
we evaluate are defined and a general strategy for one-loop calculations is summarized fol- 
lowing the method of Daniel and Sheard. ^ Results for quark bilinear operators are given 
in Sec. 3. Those for four-quark operators are described in Sec. 4 where the relation between 
lattice and continuum operators is illustrated for the case of AS = 2 operator relevant for 
extraction of the — mixing matrix. Analytical expressions for one-loop amplitudes 
are summarized in Appendices A, B and C. 

2 Formalism 

2.1 Quark operators 

The KS fermion action is given by 

S = a^Yl ^ H VM) [x{n)U^{n)x{n + fi) - x{n + /i)t/^(n)x(ri)) 

+mxin)xin) , (1) 



2a , 



where a is the lattice spacing, U^{n) denotes the gauge link variable, ri^{n) = 
and X X ^ire the single component KS fermion fields. For the construction of four- 
component Dirac fields we employ the coordinate-space method of taking a linear combi- 
nation of the X fields over a hypercube[p!2|, p!3[| . The Dirac field Q{2N) defined for each 
hypercube 2A^ G {2Z)^ is given by 

Qi2NU = WhA)^,xi2N + A) , (2) 

^ A 

where a and i are the Dirac and flavor indices, and 

7A = (3) 

with A running over the vertices of a hypercube [i.e., =0 or 1, = 1, ■ ■ ■ , 4). 
The bilinear quark operator we use is defined by 

Osf = Q{is(»^f)Q . (4) 

Here 75 = Tf^Tf^Tf^Tf** and = Ti'^^ 72 ^^73^^71'^* with the components S"^ and either 
or 1. They act on spinor and flavor indices of the Dirac field Q, and represent the spin and 
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flavor SU{4) content of tlie bilinear operator. In terms of the x fi^ld this operator can be 
written as ^ 

OsF = E Xa^s ® ^f)abXb , (5) 

AB 

where we write xa instead of x(2A^ + A) for simplicity and 

(75 ® ^f)ab = ^tr(7]i7s7B7F) • (6) 

The operator (|^) is not gauge invariant. In order to make it gauge invariant we insert 
the average of products of gauge link variables along all possible shortest paths connecting 
the cite 2N + A and 2N + B. Denoting the link factor by Uab we then define the gauge 
invariant bilinear operator by 

OsF = ^T.rAil^^)ABUfBXB , (7) 
AB 

with a and b the color indices. 

We consider two types of four-quark operators differing in the color contraction structure 
defined by 

02 = Q^ilSr®^Fr)Q''-Q\lS,®^F,)Q' , (8) 

Oi = Q''{is.®^F.)Q'-Q\is,®^F,)Q'' , (9) 



with a and b the color indices. In the method of Ref. [0] for calculating weak matrix 
elements with KS fermions the operator O2 yields an amplitude with two color contractions 
after fermion integration and Oi with a single color contraction. For this reason we shall 
refer to O2 and Oi as color two-loop and one-loop operators. These two operators generally 
mix under renormalization. 

The expression of the operators above in terms of the x fi^ld is given by 



^2 =(^y E XAils, ® ^fJabX% ■ Xcils, ® ^f,)cdXd , (10) 

^^'^^ ABCD 

^l=fl7)' E XA{lS,^^F,)ABXB-XcilS,®U)cDXD , (H) 
-■-"^ ABCD 



where we again ignored the hypercube label 2N in x(2^ + A) for simplicity. To make these 
operators gauge invariant we insert gauge hnk factors according to 

^2 =(^y E XAhs, ® ^f,)abXb ■ Xchs, ® ^F,)cDxi ■ UfnUc^D , (12) 

^-•■"^ ABCD 

= (4)' E XAilS,®^FjABXB ■ Xc{lS,®iF,)cDx'D " KdU'cB ■ (13) 
^-■-"-^ ABCD 
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For the color two-loop operators the length of paths for the link factors Uab and Ucd is 
fixed since (75- ® C,f)ab is non- vanishing only for A + B = S + F(mod2). On the other hand, 
the length varies for the factor Uad and Ucb that appear in the color one- loop operator. 



2.2 Feynman rules 

We adopt the general covariant gauge with a gauge parameter a in our perturbative calcu- 
lations. The gluon propagator is given by 

n n\ '^^A^^ n ^ •sin(afc/,/2) ■sin(aA;^/2) 

Df,^{k)ij = ^ . ^ - [1- a) ^ -2 , (14) 

and the KS fermion propagator takes the form 

f sinap^^(p -q + Tiji/a) + m5{p - q) 

J[P, -Q)ab = dab =^ : 2 ' 

where 

fi=T.^ , (16) 
J/=l 

and 

-5ip) = i2nyj2Sip+-n) . (17) 

n ^ 

The one-gluon vertex arising from the action is given by 

V^{p,-q]k)ab = -ig{T^)ab - cof^ci{p + k/2)^-5{p- q + k + TxJi/a) , (18) 

with the SU{3) generators normalized by Tr (T^T"') = 6ij/2, and the two-gluon vertex 
by 

y^il^iP, -Q; h, k2)ab = iag'^^{T^ ,T-^}ab ■ sina(p+ ^^^y^)^ 

x6^„5{p - q + ki + k2 + Tcji/a) , (19) 

where p, q are the incoming fermion momenta and k, ki, k2 the gluon momenta. 
Vertices of bilinear operators (^ have the form, 

MfkP, -<Dab = Sab^ E e^""-^-^'^'-''(¥®^)AB , (20) 
16 AB 

MskPi -Q', k)ab = -iga{T^)ab 
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(21) 



AB 



1 



X- 



16 



AB 



^iap-A—iaq-B 



{^fFW^AB{A-B)^{A-BU^^^^{aki,ak2) , (22) 



where the superscript in parentheses denote the number of emitted gluons. The function 
f(AB)i^^) defined by 



(AB) 



1 4 ,,. 

12 



(23) 



with 



1 1 



(24) 



At the one-loop level the two-gluon vertex appears only through gluon tadpole diagrams. 
Thus we only need the expression for an equal color index a = b and for the gluon momenta 
ki = —k2 = k. In this case we find 



= 1 for /i = z/ 



= e 



iak-{A^+A^) 



+ h.c for // 7^ i/ 



(25) 



with = (5 - A)^/x. 

Vertices for four-quark operators are given by a product of those for the bilinear operators 
except that the color and site indices have to be interchanged appropriately. 



2.3 Procedure of calculation 

Let us consider a one-loop diagram of a bilinear operator with two external fcrmion lines. 
The corresponding amplitude generated by the Feynman rules above are written in terms 
of momenta p taking values in the range — vr/a < p < tt/o. Since the Dirac field Q{2N) is 
defined on sites with even coordinates the physical momentum p for quarks is related to p 
through p = p + Cn/a where the vector C {Cf^ = or 1) represents the spin-fiavor content of 
quarks. We extract renormalization factors from Feynman amplitudes evaluated at vanishing 
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physical momenta p = for external fermion lines. We therefore set the external fermion 
momenta to p = Cn/a and Dtt/q. In this case the spin-flavor part of the tree amplitude 



Mg^p takes the form 



1 



(75 ® ^f)cd = - E(-l)^'''"''' ''(75 ® ^f)ab , (26) 

AB 

which shows that the calculation of renormahzation factors requires a conversion of spin- 
flavor Dirac structure from the 'single-bar' basis, in which the Feynman rules are given, to 
the 'double-bar' basis defined in (PB|). 

We employ the technique developed by Daniel and Sheard^, |^ to carry out the conver- 
sion. The general form of one-loop amplitudes that results is given by 

E / /o \d ' ^MNM'N'{ak){'jMSN ® ^M'FN')cD , (27) 



MNM'N' 



where Amnm'n' is a function of loop momenta k. The product of Dirac matrices {'Jmsn ® ^m'fn') 
can be reexpanded in terms of the basis {(75 ®^i?);S'^,F^ = 0,1}; 



ilMSN ® ^M'FN') = E Cll^M'N' ■ (7S' ® iF') • (28) 
S'F' 



The contribution of (|2^ ) to the renormahzation factor of Osf is given by 



E ^ MNM'N' / , /o \4^MNM'N'{ak) . (29) 



MNM'N' 



In most cases the decomposition (pSj) is too tedious to work out analytically. We gener- 
ate tables of Cffj^j^,^, on a computer and combined them with tables of one-loop integrals 
Amnm'n'^ separately evaluated with the Monte Carlo integration routine VEGAS, to calcu- 
late the sum (^91) . 

Our calculations are carried out for massless quark. In this case the amplitudes for the 
diagrams which have counterparts in the continuum perturbation theory contain infrared 
divergent terms of form 

d^ka^ 1 
-/a (27r)4 [E^4sin2(A;^a/2)] ■ 4 sin2(A;,a/2)] ' 

where the first factor in the denominator arises from the gluon propagator and the second 
from the massless quark propagator. To regularize the divergence we supply a finite mass n 
to the gluon propagator. The integral then takes the value 

d^ka^ 1 



(27r)4 V 4sin2(A;^a/2) ■ 4 sin2(A;,a/2) + (a/t) 
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= T7r^[-2-logK)+Foooo-7i? + l] + C'(/ta) , (31) 

with Foooo = 4.36923(1) and 7^5 = 0.577216 ■ ■ -. 

The infrared regularization above is different from that of Daniel and Sheard who added 
the mass term [anY to both the quark and gluon propagators, in which case the finite part 
of the integral is given by -foooo — 1e- We prefer not to adopt their regulator since it leads to 
a violation of fermion number conservation in continuum perturbation theory. We also note 
that the dependence on the gluon mass should cancel out between the renormalization factors 
in the continuum and on the lattice as long as one employs the same infrared regularization 
in the two cases. 

Evaluation of one-loop amplitudes for four-quark operators are much more cumbersome 
since they contain a product of two spin-flavor Dirac matrices. The calculational procedure, 
however, is essentially the same as for bilinear operators. 



3 Bilinear operators 

The one- loop renormalization of the quark bilinear operator Osf on the lattice can be written 

as 

I t 1 9^ I t 



S'F' 



where the superscript j on Ogp refer to the number of loops. The one- loop diagrams are 
shown in Fig. |I| and analytic expressions of the amplitudes are collected in Appendix A. 
Defining the coefficient z'^°^^ for the continuum operators in a similar manner the one-loop 
relation between the lattice and continuum operators are given by 



^^cont (1) _ 

'^SF — 2^ 

S'F' 



r r I y / ront lat 

OsS'OfF' + -r^{ZsF,S'F' - ZSF,S'F' 



Os'/^ . (33) 



167r^ 

The continuum renormalization factor for massless quarks can be expressed as 

Zsf!s'f' = Sss'Sff' ■ is log(/x//€) + Sss'SFF'Cr' ■ (34) 

We used the same infrared regularization as for the lattice, i.e. a finite gluon mass k, is given 
to the gluon propagator. The anomalous dimension of the operator 75 is given by 

is = 1' {CIS -I) , (35) 
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with as = (4, 1,0, 1,4) for the spin structure 75 = (/, 7^, 7^^,, 7^5, 75). The finite constant 
(jcont (jgpgnds on the continuum regularization and renormahzation schemes. For the MS 
subtraction scheme it takes the values 

^cont ^ j (10/3, 0, 2/3, 0, 10/3) for ND^ 

^ \ (14/3, 0, -2/3, 0, 14/3) for DREZ ^ ' 

for 75 = (J, 7^, 7^,^, 7^5, 75), where NDR refers to the naive dimensional regularization 
with an anti-commuting 75 and DREZ to the variant of the dimensional reduction defined 



in Ref. 15 



For the renormahzation factor on the lattice we find 

4%'F' = -Sss'Sff' ■ 75 log(aK) + C^s%s'F' ■ (37) 

The logarithmically divergent term arises from the diagrams in Fig. |l](a) and (d). It takes 
the same form for gauge invariant and non-invariant operators, and is independent of the 
gauge parameter a. Comparing (|^) and ( pT]) we find that the gluon mass k cancels out 
between Zg'p^ipi and z^gp g,p, as it should be. 

The finite coefficient C^^p gip, has the following properties. (1) The coefficients have the 
same value for the two operators with the spin-fiavor structure (75 ® ^f) and (755 ® ^p^) ( 
see Ref. |]10[ ). (2) Our explicit calculation shows that the Landau gauge part of C^gpgip, 
coming from the (sinA;'^/2 ■ sinfc'^/2) term of the gluon propagator is diagonal in spin and 
fiavor. Their values are the same for the gauge invariant and non-invariant operators. (3) 
The remaining part of the coefficient generally mixes different spin-flavor structures. Chiral 
U{1) symmetry of the KS action, however, places a restriction that operators of even distance 
with Y.1i=i{.S^ + F^){Yiiod 2) = 0,2,4 do not mix with those having odd distance (Z]^=i('S'^ + 
F^)(mod 2) = 1,3). In fact there are only a few non- vanishing off-diagonal elements (see 
Table. 0(b) below). Furthermore their values are the same for gauge invariant and non- 
invariant operators. 

Numerical values of the diagonal constants C^sFgp are tabulated in Table. |I|(a) for both 
gauge invariant (flrst column) and non-invariant (second column) operators. The results 
for the gauge non-invariant operators are given in the Landau gauge. The non-vanishing 
off-diagonal elements are tabulated in Table. |l|(b). Numerical accuracy is within 0.1%. For 
the gauge invariant operators we have numerically checked the independence of results on 
the gauge parameter a. Our results conflrm those of Daniel and Sheard ( Table 4 and 5 in 
Ref. ) after correcting our C^sfs'f' ^y —^ss'^ff'Is/'^ to take into account the difference in 
the regularization of infrared divergence. They are in a complete agreement with the recent 
calculation of the same quantity reported by Patel and Sharpe ( Table 6 and 7 in Ref. [|r^ 
). The relationship between our coefficients C^gpg,p, and theirs Cgpg,p, is 



^PS 

^SF,S'F' 



3 



6ss'Sff' (Cr* + is log tt) - C'sts'F'] , (38) 
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where C^""* are the finite continuum renormahzation factor for the DREZ scheme given by 

(H). 

We observe in Table, [^(a) that the coefficients for gauge non-invariant operators have 
a similar magnitude while those for gauge invariant ones show a substantial variation from 
operator to operator. This stems from the fact that the gauge invariant operators receive 
contribution of gluon tadpoles whose magnitude increases with the distance of the operator 
^SF = J2,i{S^ + -F^)(mod 2). Lepage and Mackenzie^ have argued that the tadpole con- 
tributions can be removed by a rescaling of fields which, for KS fermion action, takes the 
form 

X V^X , X V^X , Uq^U^ , (39) 

where Uq represents the tadpole renormahzation of link variables. A gauge-invariant choice 
for Uq is given by 



Mo 



^(Tr Up) 



1/4 1 

= l--g^ + 0{g') (40) 



12- 

with (Tr Up) the plaquette average. For gauge-invariant bilinear operators the rescaling 
amounts to a multiplication by a factor Uq~^^^. The renormahzation factors for the rescaled 
operators are obtained by subtracting 47r^(l — A5i?)/3 from the second column of Table. 0(a). 
The results listed in the third column of Table. 0(a) show that rescaled gauge invariant op- 
erators receive much less renormahzation, and that their magnitude becomes less dependent 
on the fiavor of operators. For gauge non-invariant operators without insertion of link vari- 
ables the rescaling factor is universally given by Uq. The rescaling reduces the magnitude of 
the correction without spoiling the weak fiavor dependence already apparent for the original 
operators. 



4 Four- quark operators 
4.1 Lattice result 

One-loop diagrams which contribute to the renormahzation of the color two-loop four-quark 
operators O2 defined in Sec. 2 are shown in Fig. ^ and those for color one-loop operators 
Oi in Fig. ^. In these figures horizontal lines at the four-quark vertices signify contraction 
of spin-fiavor quantum numbers, while dotted lines represent link factors and fiow of color 
indices. 

For the diagrams of Fig. |^(a)-(e) for the color two-loop operator evaluation of momentum 
and Dirac matrix parts are the same as those of the diagrams in Fig. |I| for the bilinear 
operator. The color factor is also the same for these diagrams. For the diagrams of Fig. 0(f)- 
(h), on the other hand, the color factor takes the form J2i(T^)ab(T^)a'b' which has to be 
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decomposed into the color one- and two-loop basis. This can be done by the SU (3) identity, 

'^{T^)abiT^)a'b' = -Q^ab^a'b' + ^^ab'^a'b ■ (41) 

Such a rearrangement is also generally needed for the diagrams of color one-loop operators 
in Fig. 1^ in addition to manipulation of momentum and Dirac parts. Analytic expressions 
for all the diagrams are summarized in Appendices B and C. 

Due to the mixing of color one- and two-loop operators the renormalization factor for 
the four-quark operators Oi takes a 2 x 2 matrix form, 

or = + , z,j = l,2 . (42) 

Since the operators Oi = QijSi ® (,Fi)Q ■ Q{lS2 ® Cf2)Q further depend on the pair of spin- 
flavor indices sf = (S'iFi)(S'2-F2), each element z'"* is a matrix 2;-"* = Treating 
the infrared divergence as in the case of bilinear operators we find that this matrix can be 
written as 

^ks'f = -Sff'lSss' ■ ^og{an) + C^,^^,^, . (43) 

In Table ||-^ we list the numerical values of the matrices 'jlj^sg' and C'"*^^ ,,^/ for the 
operators with the spin-flavor structure 



sf = (7M®6)(7/.®e5) , (7/.5®e5)(7/.5®e5) (Table 2 and 6) 

(7M®^5)(7/.5®e5) , (Table 3 and 7) 

(7m ® I) (7a. ® ^5) , (Table 4 and 8) 

(7^®/)(7^5®6) , (7/.5®/)(7m®6) (Table 5 and 9) 



(44) 



The results for the gauge invariant operators are in Table. and those for the gauge 
non-invariant operators are in Table The anomalous dimension matrix '^Ip^g, is gauge 
independent and takes the same value for gauge invariant and non-invariant operators. The 
results of C'-"*^^ ^/j, for the gauge non-invariant operators are for the Landau gauge. 

The one-loop renormalization coefficients for the operators having the spin-flavor struc- 
ture (755 ® ^F5)(75'5 ® ^F's) are the same as those for (75 ® C,f){is' ® C,f')- Hence the tables 
also cover the renormalization factor for the operators obtained by the interchange 7^ ^ 7^5, 
/ ^5. These operators are the most relevant for calculation of matrix elements of the ef- 
fective weak Hamiltonian. The gauge non-invariant operators with the spin-flavor structure 
(755 ® (,F5){ls' ® Cf') are renormalized in the same way as those for (7s ® ^f)(7s'' ® ^f') 
at least to one loop order, and similarly for the operators (75 ® ^f){1s'5 ® ^F'5) and (75 ® 
^f){1s' ® ^fOIIHI- known if this property persists at higher orders. We also note 

that the operators of even distance do not mix with those of odd distance due to U{1) chiral 
symmetry, similar to the case of bilinear operators. 
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The numerical accuracy is within the level of ±0.001 for the majority of elements in 
the tables, increasing to ±0.01 for large elements whose magnitude is 0(10). This accu- 
racy should be sufficient for practical applications. Reducing errors is quite computer time 
consuming because of a very large number of lattice integrals (~ 360) which have to be 
computed. 

We have checked the results in two ways, (i) For gauge invariant operators the Landau 
gauge part proportional to 1 — a has to vanish. This has been confirmed numerically, (ii) One 
can rewrite the color one-loop operator as a linear combination of color two-loop operators 
through the Fierz transformation given by 



(7s ® ^f)ab{1S' ® ^F')a'B' = — J2(^S1d ® ^E^F')AB'ilS'7D ® ^e^f)a'b ■ (45) 

DE 

The results for the color two-loop operators can then be used to evaluate the renormalization 
factors for the color one-loop operators. The results obtained in this way agree with those 
of a direct calculation of the color one-loop operators. 

The matrix elements for the mixing of gauge invariant color two-loop operators with 
color one- and two-loop operators have been computed previously by Sheard0. He took 
a basis quite different from ours, and we found it difficult to make a full comparison. For 
those cases where we can compare, however, our results are in agreement with his results. 
Also the results in Table. |] for gauge non-invariant operators agree with those of Patel and 
Sharpe[0. 

Let us finally consider improvement of four-quark operators by factoring out tadpole 
renormalizations through rescaling of fields as suggested by Lepage and Mackenzie 0. For the 
color two-loop operators the rescaling ( pQ]) yields a simple form for the improved operators, 

O-ff = (46) 

where Agf = A^^i?^ ± AS2F2 for gauge invariant operators having the spin-fiavor structure 
sf = (75i®^Fi)(752®^i^2) with AsF = ^+ Ffj,) (mod 2) the distance of bilinear operators, 
while Agf = for the gauge non-invariant operators. For the color one-loop operators, on 
the other hand, rescaling is not straightforward since the link insertion factors have lengths 
ranging from to 4. To handle this case we recall the Fierz formula ( ^5] ) and rewrite the 
color one-loop operators in terms of color two-loop operators as 

Ol;sf ='^Fsf^s>f> 02-s'f' , (47) 

s'f 

with Fsf^g/f/ numerical constants. The rescaled color one-loop operators can then be defined 



as 
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Or;'f = T.^sf,s'f'U-^''^' O,,,, . (48) 

s'f 
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For the gauge invariant choice (|40| ) for the tadpole factor uq the finite renormahzation factors 
for rescaled four-quark operators are hsted in Table. where those elements changed by 
rescaling are given after a slush (/) symbol. As one can see in the tables the rescaling indeed 
reduces the magnitude of the renormahzation correction. 

4.2 Relation with continuum operators 

In order to obtain the physical values of weak matrix elements the renormahzation factor on 
the lattice obtained in the previous section has to be combined with those in the continuum. 
In this section we illustrate the procedure for the K meson B parameter relevant for the 
— mixing matrix. 

The K meson B parameter Bk in the continuum theory is defined by 

= 02^-2 • (49) 

3J K"''K 

In the method of Ref. Ijl^ for calculating weak matrix elements with KS fermions, the 
operator in the numerator is replaced by the sum of the following four operators 

A2 = S'^ii.^ ® ^5)/^" ■ S\-f,5 ® ^5)D' , (50) 

where S and D are the KS quark fields introduced for s and d quark separately, a and b the 
color indices, and the quark fields in the first current are to be contracted with and those 
in the second current with K^. The choice of fiavor ^5 in these operators corresponds to the 
use of -D(75 ® ^5)5' for creating the external and in the Nambu-Goldstone channel 
associated with U{1) chiral symmetry of the KS fermion action. 

The renormahzation factor for these operators can be read off from Table. 2 for gauge 
invariant operators and from Table. 6 for gauge non-invariant operators. As can be seen, the 
four operators not only mix among themselves but also with a large number of others having 
different spin-flavor structures. We note that the extra operators all have the flavor matrix 

7^ ^5- Since the and mesons are created with the flavor ^5, the matrix element 
of the extra operators should vanish in the continuum limit where a restoration of SU{4) 
flavor symmetry is expected. The renormahzation factor for some of the extra operators are 
numerically not small, however. Whether they yield negligibly small contributions at the 
current range of inverse lattice spacing 1/a ~ 2 — 3GeV has to be checked through actual 
simulations. For simplicity we disregard the mixing with the extra operators in the following. 
Extensions to the general case is straightforward. 
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Denoting the four operators ( |50D as {O 



lat . 



a 



4x4 anomalous dimension matrix 7^^^ is given by 



lat 

laf3 



( 9 


-7 
V -6 



-3 


-3 
2 



-7 
-6 
9 




, 4} = {Vi, ■ ■ ■ , ^2}, we find that the 



-3\ 
2 

-3 

/ 



(51) 



which takes the same form for the gauge invariant and non-invariant cases. The finite part 
takes the values, 



-18.915 


-4.772 


-5.253 


-2.251 \ 





-60.000 


-4.502 


1.501 


-5.253 


-2.251 


-19.513 


-2.977 


-4.502 


1.501 





/ 



(52) 



for the gauge invariant operator, and 



( 37.446 


-5.253 
\ -4.502 



-2.913 
28.706 
-2.251 
1.501 



-5.253 
-4.502 
37.976 




-2.251 \ 

1.501 
-4.504 
24.464 / 



(53) 



for the gauge non-invariant operator in the Landau gauge 
invariant operator (|5^) is in agreement with those of Ref. ||lT 
row of 



The result for the gauge non- 
The second and the fourth 



_)2D has previously been calculated by Sheard, with which our results agree. For the 
rescaled operators discussed in Sec. 4.1 the diagonal elements of C^'^ changes to ( 7.403, 
—7.361, 6.805, ) for the gauge invariant operators and to ( 11.127, 2.387, 11.657, —1.854 ) 
for the gauge non-invariant operators (off-diagonal elements are not affected). 
The IS.S = 2 continuum operators are given by 



C2 



with the one-loop renormalization taking the form 

,2 



(1) 



AO) 
^3 



(54) 



(55) 



where the anomalous dimension matrix '^1° is given by 



cont 
hj 




(56) 
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For the finite part Cf?" we find that 



f^cont cont 



with 



11/12 
7/12 



for NDR[17] 
for DREZ 



(57) 



(58) 



In order to relate the lattice operators to those in the continuum let us define a 2 x 4 
matrix by 

1 1 0~ 
10 1 



M 



(59) 



with which one can write C 



(0) 



MiaO'^^ Using the intertwining property Mi^'y^^p 



Mjp it is easy to see that the one-loop renormalization relation between the continuum 
and lattice operators is given by 



'(1) 



M,, 



Sap + ^(V"* log(/.a) + c^^^' - C^^X^ 



O 



lot (1) 



(60) 



The numerator of the parameter with one-loop renormalization correction equals the 
sum (i?0|4'^|i^°) + (i?°|4'^|i^°)- 
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Appendix A One-loop amplitudes for bilinear operators 

The one-loop amplitudes corresponding to the diagrams of Fig. |l] for the external quark 
momenta p = Cn/a, Dir/a and color indices a, b are as follows(the common factor 4/3 ■ 6ab ■ 
(^^/(IGvr^) is not included). 



{as - (1 - «)) a;(75 (g) ^f)cd + ^mn {it^pMSNa^, ® ^mfn)cd 

^ipaMN 

(1 - a) ^ Xm{imsm ® ^mfm)cd 

M 

^MnI^ + ^\^J.{l^lp^l5MSN ® C,h5MFN + 1^51^1 SNp^i ® C,p5Mfn)cd 

MNfip 

MNp 



-(3 + a)Zoooo(7s ® ^f)cd 
-xa — -(1 + q;)Zoooo — -R 



(75 ® i 



F)CD 



--(3 + a)A5j.Zoooo - 2(1 - «)TAg^ 



(7s ® ^f)cd 

21ogaK + Foooo-7£;+l, l'S' + F|^ 



(61) 

(62) 

(63) 
(64) 
(65) 



where as = (4, 1, 0, 1, 4) for S = (J, 7^, 7^^, 7^5, 75), x = 
(S^ + F^) (mod 2), and A^f = E;. + i^U- 

We use the following notations to simplify the expressions for the loop integrals which 
appear in the amplitudes above: 



= sm (pf, , 



= sin 0^/2 , B=[4j2sl^'' 



Cf, = cos (p^ , 



COS (j)f,/2 



,21-1 



167r^ 



(27r)^ 



(66) 



where is the loop momentum and B and F originate from gluon and fermion propagators. 
In terms of these symbols the loop integrals are defined as follows. 



X^N = I \clspS,EM{<P)EN{-(t>)BF^ - UpJu^^mB^' 
J (h L 4 



X 



Em{(P)Em{-<I)) - 6 



(67) 
(68) 

(69) 



j = l o-^/i 
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j = l (Jy^^ 



Zqqoo —IB 

Jcf, 



Ta= f 2slslB^ X (0, 0,1,2 + C3, 3 + 2c3 + C3C4) 
j(f, 

for A = (0,1,2,3,4) 
R = I \c,{2sl - i)(-2 - 2sl + ^)BF' - B' 



where and O^^i, are given by 



Em{<P) =l[l + i-lf-e'^-/') 



(70) 
(71) 

(72) 
(73) 



(74) 
(75) 

(76) 



p=i 



p=i 



The integrals X^^f^ and iX'^v "were evaluated numerically by Daniel and Sheard]^. We 
have confirmed their results except for some sign reversals for Yj^j^ in their Table. 3. For 
the calculation of the integral we employed the Monte Carlo integration routine VEGAS. 



Appendix B One-loop amplitudes for color two-loop four- 
quark operators 

Analytic expressions for one-loop diagrams in Fig. ^ are listed below for color two-loop four- 
quark operators of a general spin-flavor structure (7s®Cf)(75' ®^f')- The amplitudes for the 
diagrams Fig. |^(a)-(e) are products of tree and one-loop bilinear amplitudes. The diagrams 
in Fig. |](f)-(h) cannot be factorized in this way. They take the following form, where we 
drop the common factor J2i{T^)ab{T^)a'b' and ^f^/lGvr^. External fermion lines have momenta 
p = Cn/a-, Dn/a, Cn/a, D'n/a and color indices a,b,a',b' as specified in Fig. ^. 
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g2(^)= 

MP 



~ X! ^Mn{1ppMS ® iuF — ISMpp® iFM)cD 
fipaMN 



X { ipcrMS' ® ^M F' — IS'Map ® ^F' m)c'D' 

^ ^m{1MS ® ^MF — Ism ® iFM)cD{,lMS' ® C,MF' — IS'M ® C,F'Al)c'D' 
M 

(77) 



fj.pLMN 



D' 



\S + -^lAt(7/i5MSL ® ip^MFLjCoilS'Npp ® ^F'N — IfipNS' ® ^NF')c' 
+ \S' + F'\^{-fSNpp®^FN 



[l-a] 

pLMN 



+ \S' + F'\^{jsN ® 'CfA' - 7AfS ® ^Nf)cd{1p5MS'L ® ipbMF'LjC'D' 



(7J 



~ X! ^^LMAfl'S' + -^UI'S" + -F'|^(7At5-ft:5L ® ^p5KFl)cd{1 p5MS' N ® ^p5MF'n)c'D' 
fiKLMN 

+ (1 - ") H ^XLAfTvl'^ + i^U(7/.5i^SL ® ^p5KFl)cD 
puKLMN 

X\S' + F'\y{jy^MS'N ® ^u5MF'n)c'D' 

(79) 

where jS" + F\^ = {S^ + F^)(mod 2). The four types of loop integrals are defined by 

U'lU = J ^c,s,BF^ Y: E El{91^:^)Em{4> - e^})EM{-<P) (80) 
U^,,^ = f i2-s,B'^ E E,{e^})EM{<P - ef^)E^{-ct^) (81) 

•'^ v+pj=\ 

vLmn = j B^Y.ilEK{e^E:^{<p-e^ 
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X 1^ E E EM{-elibEN{-<p + ef,) (82) 

4 



Xl^EE^MM^^)i?^(-0 + ^l^) . (83) 

The numerical values of the last three integrals reported in Tables 4 and 5 in Ref. have 
some minus signs missing. We have corrected the sign and evaluated the values of Uj^fj^ 
which are not listed in Ref. 0. 

Appendix C One-loop amplitudes for color one-loop four- 
quark operators 

Analytic expressions for one-loop diagrams in Fig. |^ are listed below for color one-loop four- 
quark operators of a general spin structure (75 ® C,f){1s' ® ^f')- We set the momenta of 
the external fermion lines to p = Cn/a, Dn/a, Cn/a, D'n/a and assign the color indices 
a, b, a', b' as shown in Fig. |[ The amplitudes of the diagrams Fig. |5(a), (c), and (d) are the 
same as those for the color two-loop operators except for the color factor which takes the 
form Y^j(T^)ab'(T^)a'b for the diagram (a) and Sab'Sa'b for the diagrams (c) and (d). Other 
amplitudes take the following form, where we drop the common factor g'^/lGir'^. 



{lupMS ® iMF)cD{lS'N ® ^F'n)c'D' + {iSMpp ® ^Fm)cd{'JNS' ® (.NF')c'D 
+ {lSN ® ^FN)cD{lfMpMS' ® ^MF')c'D' + {InS ® f) C d{.1 S' M pp ® ^F'm)c'D' 

-4(1 - «) E>^£5M,M ((-l)^^-(^--) + (-l)--(^'+-')) 
pM 



X (75M ® ^fm)cd{7s'm ® ^f'm] 



CD' 



(84) 



-2^0000 (7s ® ^f) coils' ® i 



F'jC'D' 
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+ (1 - a)-Zoooo{ls ® ^f)cd(7s' ® ^f')c'D' 

-(1 - a)5]Zoooo^((-l)^''+^^'* + (-l)(^'+^')-)(7^55 ® W)cd(7m55' ® W)c'D' 



_((_l)Af-(5+F) ^ (_i)M-(5'+F'))(^^^^ ^ eMF)cD(7M5' ® ^MF'h'D' 

x(7a»5MS ® i|J.bMF)cD{l^lbMS' ® inbMF')c'D' 
_((_1)(-S'+-^)m+(-S'+-F)^+M-(5+F) _^ j^_]_^(5'+F')m+(5"+F')^+M-(5'+F')^ 



xilfMuMS ® (.fiuMF)cD{lfiuA4S' ® (.iiuMF')c' 



D' 



UpaMN 



illipMS ® iMF)cD{lS'Ncjp ® ^F'n)c'D' + ilSMpp ® ^Fm)cd{1 pcjN S' ® ^NF')c> 

I 



ilfipMS ® C.MF)cD{lpaNS' ® ^NF')c'D' + {iSMpp ® ^Fm)cd{1S' Nap ® 'CF'Ar)c" 



D' 



D' 



+ (1 - E(^A^ + ^'^Mo)(75M ® ^FM)cDilS'M ® ^F'm)c'D' 



M 



15) 



(86) 



E r)i^[fl5L]MN ^[p.5M]LN> 
ppLMN ^ 



{iLSNpp ® (.lfn)cd{1s'm ® (.f'm)c'd' + (75M ® (,fm)cd{1ls'Npp ® ^LF'n)c'D' 

+ {lppNSL ® ^NFL)cDilMS' ® ^MF')c'D' + (ImS ® ^MF)cDilppNS' L ® ^NF'l)c'D' 
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liLMN 



2 \'-^[fi5L]MN '^[^ibM]LN 



{iSM ® ^fm)cd{1ls'n ® ^lf'n)c'd' + ilMS ® ^mf)cd{1ns'l ® (,NF'l)c'D' 

+ {lLSN ® ^LFn)cd{1S'M ® C,F'm)c'D' + {iNSL ® C,NFl)cd{1MS' ® C,MF')c"D' 



T(7i<:sAf ® C,kfn)cd{1ms'l ® C,mf'l)c'd' 



{3+F+3'+F')^ 



,,,, li' r »,f ^ ^ 



)^/^5-ff]i[At5M]Af 



The new integral Tj} is defined by 



where 



the components /i, i/, p, and a are all different. 
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The operator E is independent of the left-left operator F ® F in general dimensions, 
vanishing in four dimensions. It can be ignored at one loop, but has to be retained at 
two loop order and beyond. 
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Tabic 1: One- loop perturbative corrections for bilinear operators. Operators of form (75 <S) 
Cf) and (755 ® C-Fs) receive the same one-loop corrections. The components fi, p and a 
are all different and not summed. Values of anomalous dimension 75 are also listed, (a) 
Diagonal elements of the renormalization factor for gauge invariant (first column) and non- 
invariant (second column) operators. The third and fourth column shows the values for 
rescaled operators as discussed in text, (b) Off-diagonal elements of the correction which 
take the same values for gauge invariant and non-invariant operators. The first and second 
columns of operators specifies the row and column of the mixing matrix. 



(a) 


operator 


7 


invariant 


non-invariant 


inv (rescaled) 


non-inv (rescaled) 


1 


(/® 


I) 


8 


55.585 


55.585 


42.426 


42.426 


2 


[1® 


C5) 


8 


-47.783 


12.813 


-8.304 


-0.346 


3 


[l® 


Cm) 


8 


14.844 


27.077 


14.844 


13.918 


4 


[l® 


Cms) 


8 


-29.948 


14.405 


-3.629 


1.246 


5 


[l® 


Cmi') 


8 


-10.569 


17.583 


2.589 


4.423 


6 


(7m ^ 


3/) 





0.000 


12.232 


0.000 


-0.927 


7 


(7m ^ 


5Cs) 





-30.000 


14.353 


-3.682 


1.193 


8 


(7m ^ 


5 Cm) 





19.693 


19.693 


6.533 


6.533 


9 


(7m <^ 


5C.) 





-13.388 


14.764 


-0.228 


1.605 


10 


vO'm 


5 U) 





-13.409 


14.743 


-0.249 


1.584 


11 


(7m ® Cms) 





-45.988 


14.608 


-0.651 


1.448 


12 


(7m ^ 


3 Cm-) 





4.519 


16.752 


4.519 


3.593 


13 


(7m ^ 


5C-a) 





-29.651 


14.702 


-3.332 


1.543 


14 


v7m'^ 


®I) 


-8/3 


-14.623 


13.529 


-1.464 


0.369 


15 




®Cm) 


-8/3 


-0.428 


11.804 


-0.428 


-1.355 


16 




®Ca) 


-8/3 


-29.668 


14.685 


-3.349 


1.525 


17 




®Cm-) 


-8/3 


7.728 


7.728 


-5.430 


-5.430 


18 




^Cma) 


-8/3 


-14.200 


13.952 


-1.041 


0.792 
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®^Xa) 


-8/3 


-45.390 


15.206 


-5.911 


2.046 



(b) 

operator mixed operator 



9 


(7m € 




(7m 


5 Cm) 




-4.504 


11 


(7m 


5 Cms) 


(7m S 


5C.5) 




0.860 


13 


(7m 




(7m ^ 


5 Cm-)' (7m 


^Cam) 


1.980 


16 


(7m- 


»Ca) 


(7m- 


»Cm)' (7m- 


8C-) 


0.902 
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Table 2: Anomalous dimensions and finite corrections for gauge invariant four-quark oper- 
ators with the spin-flavor structure VV = (7^ ® 'C5)(7m ® Cs) = (7^5 ® ^5){lii5 ® ^5)- 
Subscripts 1 or 2 attached to the operators refer to the the number of color loops. Operators 
listed at the top row mixes with those on the first column with the numerical coefficients 
given (a common factor g'^/ldn'^ is removed). The values after '/' are the finite corrections 
for the rescaled operators. The two rows of numerical values for each operator in the first 
column are for the color one-loop operator (first row) and for the color two-loop operator 
(second column). All indices of operators are summed with different indices not taking equal 
values. 





mixed operator 




7 


finite 


7 


VV2 
finite 


7 


AAi 
finite 


7 


AA2 
finite 


1 


(Z® 








4.725 








-19.114/ - 5.954 




-5.573 












-1.575 








-2.412 




1.857 


2 


(-f ® 


eM){J®e-) 






-0.624 








-0.485 




1.112 












0.208 








-0.208 




-0.370 


3 


(-f ® 


ea6)(-f ®ep5) 




-0.052 




-0.235 




-0.067 
















-0.022 




0.078 




0.022 






4 


(-f ® 


ea6)(-f ^Cas) 




-0.876 




-1.310 




0.565 
















-1.029 




0.436 




-0.188 






5 


(75 i 


?>Cm)(75 ®e 


m) 




20.774/7.615 








-3.675 




-3.150 












-2.569 








-1.575 




1.050 


6 


{75 (5 


Km){75 




0.319 








0.485 




0.416 












0.706 








0.208 




-0.138 


7 


{75 (5 






0.052 




0.044 




0.130 
















0.022 




-0.014 




-0.213 






8 


{75 ^ C<t5)(75 <S> 






-0.439 




-0.376 




1.511 
















-0.188 




0.125 




-0.872 






9 


(7m i 


S>/)(7m ® ^) 






0.428 








-0.086 




-1.107 












-0.142 








-0.396 




0.369 


10 


{7m i 


s>6){7m ®6) 


9 


-18.915/7.403 






-7 


-5.253 


-6 


-4.502 










-3 


-4.772 




-60.000/ - 7.361 


-3 


-2.251 


2 


1.501 


11 


{7m i 


?>?mi')(7m «> 






-1.289 




-1.104 




2.683 
















-0.552 




0.368 




-1.615 






12 


{7m i 


S>$mi')(7m «> 


?Mp) 




0.174 




0.149 




0.390 
















0.074 




-0.049 




-0.573 






13 


(7m « 


3 4m"')(7m ® 






-0.051 




0.402 




















0.153 




-0.134 










14 


(7m < 


?l?i'lt,)(7M ® 






-1.139 




-2.774 




2.342 
















-1.189 




0.924 




-0.780 






15 


(7m « 


?i^./|p)(7m ® 






0.277 




0.149 




0.224 
















-0.232 




-0.049 




-0.074 






16 


(70-5 


®-'")(7o-5 ® 


I) 




-0.086 




-0.478 




0.428 
















-0.396 




0.159 




-0.142 






17 


(7ctB 


® ^b)(7ct6 @ 




-7 


-5.253 


-6 


-4.502 


9 


-19.513/6.805 














-3 


-2.251 


2 


1.501 


-3 


-2.977 






18 


(7ctB 


® ^Mi/)(7ct5 






2.342 








-1.007 




-2.774 












-0.780 








-1.584 




0.924 


19 


(7ctB 


® ?M'')(TctB 






0.224 








0.174 




-0.656 












-0.074 








0.074 




0.218 


20 


(70-5 


® ^M"')(Tct5 






0.083 
























-0.249 














21 


(7ctB 


® ^M<^){7ct6 


®?M<^) 




2.683 








-1.289 




-1.104 












-1.615 








-0.552 




0.368 


22 


(7ctB 


® ^Mct)(7ct6 






0.224 








0.174 




0.149 












-0.074 








0.074 




-0.049 


23 


(7m^ 


® ^m)(7m"' i 


3«m) 




-3.224 




-2.763 




21.486/8.327 
















-1.381 




0.921 




-3.157 
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(7m"' 


i8i?m){7m"' « 






-0.485 




-0.416 




-1.465 
















-0.208 




0.138 




2.733 






25 


(7mi' 


(X>5m)(7mi' ^ 


9€p) 




-0.222 




0.764 




















0.668 




-0.254 











25 



26 




-19.434/ - 6.274 


-5.960 


4.145 








-3.000 


1.986 


-1.381 




27 




0.485 


1.944 


0.624 








0.208 


-0.648 


-0.208 




28 


® ?i/5){7mi' ® Cms) 


-0.003 




-0.052 


-0.044 






-0.168 




-0.022 


0.014 


29 


(7Mf ® ?i'5)(7mi' ® Ci^s) 


0.918 




-0.900 


-0.771 






-0.675 




-0.385 


0.257 


30 




0.031 












-0.095 








31 


(7mi' 'X' 5<t5)(7mi' ® 5p5) 


0.067 




0.052 


-0.528 






-0.022 




0.022 


0.176 


32 


(7^^ ® €a5)(7Mi' ® Cas) 


1.157 




-0.416 


-1.704 






-0.385 




-0.832 


0.568 
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Table 3: Renormalization factors for the gauge invariant operators VA = (7^(8>^5)(7^5®.^5). 
Matrix elements are arranged in the same way as in Table. ^ 





mixed operator 




7 


V ^1 
finite 


7 finite 


T 
1 




) ^)(7m5 8) 


^) 


















—0.142 




9 


(7n « 


)€5)(7m5 € 




Q 


— ly.oiy/D.yyo 












—3 


—3.560 


— 29.999/ — o.DoU 


Q 
O 


(7m « 


)?M>')(7m5 


i^Cm-') 




— i.zoy 














—0.552 


0.368 


A 


(7m '8)?mi')(7m5 iXi^Mp) 






u. i^y 












0.074 


—0.049 


r: 
O 


(7m S 


)Cm>')(7m5 


® ?^p) 




— U.UOl 














0.153 




U 


(7m S 


)f^p)(7M5 


fXl^Mf) 




















—0.134 


7 
1 


(7m ® ?i'p)(7m5 


fXiCi-p) 






1 1 ^9 












—0.576 


0.384 


8 


(7m S 


)f^p)(7M5 


(X>?i.<t) 




0.225 


-0.253 












-0.078 


0.084 


9 


(7ct5 ' 


X>-f)(7<T ® 


I) 




-0.191 


-0.164 












-0.082 


0.054 


10 


ilrrS ' 


X>6)(7<T 


5 6) 


-7 


-5.253 


-6 -4.502 










-3 


-2.251 


2 1.501 


11 


(7<t5 < 


X>Cmi')(7<t 


«'?M'') 




2.342 














-0.780 




12 


ilrrS ' 


X>Cmi')(7<t 


■^^MP) 




0.224 














-0.074 




13 


(ItS ' 


X>Cmi')(7<t 


«'?M't) 




0.083 














-0.249 




14 


(7ct5 ' 


Cmo-)(7o- ® Cmct) 




2.412 














-0.804 




15 


(7<t5 ' 


X>Cm<t)(7<t 


® 6<t) 




0.307 














-0.324 




16 


(7MI' 


»5m)(7p<t 


® 65) 




0.174 


0.149 












0.074 


-0.049 


17 


(7m>' 


X)Cp)(7p<t 


® 65) 




-0.091 


0.253 












-0.324 


-0.084 


18 


(7mi' 


^ 5i'5)(7p<T Cp) 




0.173 














0.078 




19 


(7m>' 


» Cct5)(7p<t ® Cm) 




-0.224 














0.074 
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Table 4: Renormalization factors for the gauge invariant operator VV = (7^ ® I){'Jfi ® ^5)- 
Matrix elements are arranged in the same way as in Table. |[ 

VVi VV2 



mixed operator 7 finite 7 finite 



1 


{■y^, (S) I)i■y^, <S> is) 


9 


— 19.319/6.999 










—3 


—3.560 


— zy.yyy/ — o.bou 


2 


{7m«>?5){7m'X'/) 




0.428 












n 1 /to 




3 


{7m 'S'5mi^)(7m ® €i/<t5) 




0.051 












-0.153 




4 


{7m 'S)5m^)(7m ® €m<t5) 




-0.225 


0.253 










U.U / 


— U.Uo4 


5 


{7m « 






— 1.344 


— 1.152 










— U.O/D 


U.oo4 


6 


{7m S 


)?^p)(7m <»S.p^5) 




0.174 


0.149 










U.U / 4 


— U.U4y 


7 


{7m « 


)C^p)(7m 'XjCi.ps) 




— 1.289 


— 1.104 










— U.ooz 


U.oDO 


Q 



{7m « 


)Ci.p)(7M iXl^Mps) 


















-0.134 


9 


{7<t5 1 




-7 


-5.253 


-6 -4.502 








-3 


-2.251 


2 1.501 


10 


{7<t5 < 


X>C5){7<t5 




-0.191 


-0.164 










-0.082 


0.054 


11 


{7<t5 < 


Xl Cpi.){7<T5 ® (jyab) 




-0.083 












0.249 




12 


{7<t5 < 


2> 5mi'){7<t5 ® Ci'ps) 




0.307 












-0.324 




13 


{7<t5 < 


X> Cmi'){7<t5 ® Cmi-s) 




2.412 












-0.804 




14 


{7<t5 ' 


8?M<^)(7<t5 ^^pcrs) 




-0.224 












0.074 




15 


{7<t5 I 






2.342 












-0.780 




16 


{7<t5 < 


X'Cpct)(7<t5 (^Cmps) 








17 


{7MI' 


8>Cm)(7p<t5 8>5p) 




0.091 


-0.253 










0.324 


0.084 


18 


{7m"^ 


^Cp)(7p<t5 (^^m) 




-0.174 


-0.149 










-0.074 


0.049 


19 


{7m"' 


^?i'5){7p<t5 IXI^cts) 




0.224 












-0.074 




20 


{7mi' 


8>Cct5){7p<t5 fXl^i^s) 




-0.173 












-0.078 
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Table 5: Renormalization factors for the gauge invariant operator VA = (7^ (g) I){pifj.^ ® ^5) 
and AV = {"^^^ ® ® ^5)- Matrix elements are arranged in the same way as in Table. 0. 













VA2 




AVi 


AV2 




mixed operator 


7 


finite 


7 finite 


7 


finite 


7 finite 


1 




?m)(75 ^i^ms) 




1.511 
-0.872 






-0.439 
-0.188 


-0.376 
0.125 


2 




eM)(75 «>Ca5) 




0.130 
-0.213 






0.052 
0.022 


0.044 
-0.014 


3 




C<t5){75 ®Cm) 




0.485 
0.208 


0.416 
-0.138 




0.319 
0.706 




4 


(I» 


C<t5){75 




-3.675 
-1.575 


-3.150 
1.050 




20.774/7.615 
-2.569 




5 


(75 i 


5CM){/®eM5) 




0.565 
-0.188 






-0.876 
-1.029 


-1.310 
0.436 


6 


{75 i 


5CM){/^?a5) 




-0.067 
0.022 






-0.052 
-0.022 


-0.235 
0.078 


7 


(75 i 






-0.485 
-0.208 


1.112 
-0.370 




-0.624 
0.208 




8 


(75 i 


5C<t5)(/®C<t) 




-19.114/ - 5.954 
-2.412 


-5.573 
1.857 




4.725 
-1.575 




9 


(7m <: 


g>/)(7M5 ® 6) 


9 


-19.513/6.805 




-7 


-5.253 


-6 -4.502 








-3 


-2.977 




-3 


-2.251 


2 1.501 


10 


{7m <: 


g>6){7M5 f^-f) 




0.428 
-0.142 






-0.086 
-0.396 


-0.478 
0.159 


11 


{7m i 


g'CMi')(7M5 (^?i.<t5) 










0.083 
-0.249 




12 


{7m <: 


Cmi')(7m5 iX) Cmcts) 




-0.174 
-0.074 


0.656 
-0.218 




-0.224 
0.074 




13 


{7m <: 


Cmi')(7m5 (Xl 




-1.007 
-1.584 


-2.774 
0.924 




2.342 
-0.780 




14 


{7m i 


g) Ci^p){7m5 ® Spas) 




0.174 
0.074 


0.149 
-0.049 




0.224 
-0.074 




15 


{7m <: 


Ci^p){7m5 ® Cfps) 




-1.289 
-0.552 


-1.104 
0.368 




2.683 
-1.615 




16 


{7m i 


?'Ci^p){7m5 ® Cmps) 










-0.083 
0.249 




17 


{7<t5 


iX>/)(7ct ® 6) 


-7 


-5.253 


-6 -4.502 


9 


-18.915/7.403 










-3 


-2.251 


2 1.501 


-3 


-4.772 


-60.000/ 


18 


{7<t5 


«>C5)(7^ 




-0.086 
-0.396 


-1.107 
0.369 




0.428 
-0.142 




19 


{7<t5 


<X> Cmi'){7<t ® Si^as) 










0.051 
-0.153 


-0.402 
0.134 


20 


{7<t5 


<X>Cm"^){7<t ®?i^p5) 




0.390 
-0.573 






0.174 
0.074 


0.149 
-0.049 


21 


{7<t5 


fX Cmi'){7<t ® Smi^s) 




2.683 
-1.615 






-1.289 
-0.552 


-1.104 
0.368 


22 


{7<t5 


lX> Cmct)(7ct (X Cpcrs) 




-0.224 
0.074 






-0.277 
0.232 


-0.149 
0.049 


23 


{7<t5 


lX> Cmct)(7it iX Cmo-s) 




2.342 
-0.780 






-1.139 
-1.189 


-2.774 
0.924 


24 


{7<t5 


(X>Cm<^)(7<t fXffips) 










-0.051 
0.153 


0.402 
-0.134 


25 


{7mi' 


iX Cm)(7mi^5 <8) Cms) 




-0.416 
-0.832 


-1.704 
0.568 




1.157 
-0.385 





29 



26 




0.052 


-0.528 


0.067 








0.022 


0.176 


-0.022 




27 








0.031 












-0.095 




28 








0.031 












-0.095 




29 




-0.900 


-0.771 


0.918 








-0.385 


0.257 


-0.675 




30 




-0.052 


-0.044 


-0.003 








-0.022 


0.014 


-0.168 




31 




0.624 




0.485 


1.944 






-0.208 




0.208 


-0.648 


32 




4.145 




-19.434/ - 6.274 


-5.960 






-1.381 




-3.000 


1.986 


33 








-0.222 


0.764 










0.668 


-0.254 


34 


ilniy ® ?CTB)(7mi/5 ® Cm) 






-0.222 


0.764 










0.668 


-0.254 


35 


(7mi' ® CCTB)(7fii^5 ® Cp) 


-1.465 




-0.485 


-0.416 






2.733 




-0.208 


0.138 


36 


(7mi' ®€a5)(7Mi^5 ® Ccr) 


21.486/8.327 




-3.224 


-2.763 






-3.157 




-1.381 


0.921 
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Table 6: Renormalization factors for the gauge non-invariant operators VV = (7^ (8)^5) (7^® 





and A A = (7.5 (g) ^5) (7.5 


(g) ^5) in the Landau gauge. 


Matrix elements 


are arranged m 


the 


same way as in Table. | 


1 


















VVi 


VV2 


AAi 


AA2 




mixed operator 


7 


finite 


7 finite 


7 finite 


7 finite 


1 




Cm)(^®«m) 




4.725 




-3.774 


-3.235 










-1.575 




— 1.617 


1.078 


2 


(-fig 


Cm)(/ (g?!.) 




-0.624 




-0.485 


-0.416 










0.208 




—0.208 


0.138 


3 


{-f (g) 






-0.052 


-0.044 


-0.067 












-0.022 


0.014 


0.022 




4 


(/ Ig 


C<T5)(/<gC<T5) 




-1.466 


-1.257 


0.565 












-0.628 


0.419 


—0.188 




5 


{75 (5 


5Cm)(75 fg^M) 




4.853 




-3.675 


-3.150 










-1.617 




— 1.575 


1.050 


6 


(75 5 


??m)(75 




0.624 




0.485 


0.416 










-0.208 




0.208 


—0.138 


7 


(75 5 


K<t5)(75 (gSps) 




0.052 


0.044 


0.067 












0.022 


-0.014 


—0.022 




8 


{75 (5 


5C<t5)(75 (g^CTs) 




-0.439 


-0.376 


1.885 












-0.188 


0.125 


—0.628 




9 


(7m C 


g/)(7f. (g/) 




0.428 




-0.538 


-0.461 










-0.142 




—0.230 


0.153 


10 


(7m i 


g6)(7M 8 6) 


9 


37.446/11.127 




-7 -5.253 


-6 -4.502 








-3 


-2.913 


28.706/2.387 


—3 —2.251 


2 1.501 


11 


(7m <: 


gCMi')(7M 'g^Mi') 




-1.289 


-1.104 


3.016 












-0.552 


0.368 


— 1.005 




12 


(7m <: 


gCMi')(7M (g^Mp) 




0.174 


0.149 


0.224 












0.074 


-0.049 


—0.074 




13 


(7m i 


gCi^p)(7M gl^i^p) 




-1.813 


-1.554 


2.342 












-0.777 


0.518 


—0.780 




14 


(7m <: 


g?>^p)(7M 'g^i^.r) 




0.174 


0.149 


0.224 












0.074 


-0.049 


—0.074 




15 


(7<t5 


'g^)(7<T5 fg-f) 




-0.538 


-0.461 


0.428 












-0.230 


0.153 


—0. 142 




16 


(7<t5 


'gC5)(7<T5 (gCs) 


-7 


-5.253 


-6 -4.502 


9 37.976/11.657 










-3 


-2.251 


2 1.501 


—3 —4.504 


24.464/ — 1.854 


17 


(7<t5 


<g Cm"^)(7<t5 ® Cm'') 




2.342 




-1.813 


-1.554 










-0.780 




—0.777 


0.518 


18 


(7<t5 


<g Cmi')(7<t5 gl Cmp) 




0.224 




0.174 


0.149 










-0.074 




0.074 


—0.049 


19 


(7<t5 


Ig Cmct)(7ct5 <g Cmct) 




3.016 




— 1.289 


— 1.104 










-1.005 




—0.552 


0.368 


20 


(7<t5 


Ig Cmct)(7ct5 Ig Ci'o-) 




0.224 




0.174 


0.149 










-0.074 




0.074 


—0.049 


21 


(7m<^ 


lgCM)(7M>' -g^M) 




-3.224 


-2.763 


5.433 












-1.381 


0.921 


— 1.811 




22 


(7m<^ 


lgCM)(7M>' -gC-^) 




-0.485 


-0.416 


—0.624 












-0.208 


0.138 


0.208 




zo 


(7mi' 


lgCp)(7Mi' ®5p) 




-4.226 


-3.622 














-1.811 


1.207 


-1.381 




24 


(7mi' 


'gep)(7Mi' ®C<t) 




0.485 


0.416 


0.624 












0.208 


-0.138 


-0.208 




25 


(7MI' 


lg?^5)(7Mi' (g^Ms) 




-0.067 




-0.052 


-0.044 










0.022 




-0.022 


0.014 


26 


(7MI' 


Ig Cj/5)(7mi' Ig Ci^s) 




1.293 




-0.900 


-0.771 










-0.431 




-0.385 


0.257 


27 


(7MI' 


<g C<t5)(7mi' ® Cps) 




0.067 




0.052 


0.044 










-0.022 




0.022 


-0.014 


28 


(7m<^ 


Ig C(t5)(7mi' ® Co-s) 




1.157 




-1.005 


-0.862 










-0.385 




-0.431 


0.287 
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)le 7: Renormalization factors for the gauge non-invariant operators VA = (7^®^5) (7^ 
in the Landau gauge. Matrix elements are arranged in the same way as in Table. 0. 



VAi VA2 



mixed operator 7 finite 7 finite 



1 


{7m « 


5^)(7m5 




0.428 
-0.142 





z 


{7m ^ 


S ?5j(,7M5 ?5j 


y 

-3 


QV 71 1 /II QOO 

( . ( 11/ ll.vSyz 

-3.708 


26.585/0.266 


3 




5Cm"^){7m5 ® Cm"^) 




-1.289 
-0.552 


-1.104 
0.368 


4 


{7m « 


5Cmi'){7m5 i^Cmp) 




0.174 
0.074 


0.149 
-0.049 


5 


{7m « 


5Ci^p){7m5 ® Ci-p) 




-1.813 
-0.777 


-1.554 
0.518 


6 


{7m « 


5 Cixp){7m5 ® ^l-cr) 




0.174 
0.074 


0.149 
-0.049 


7 


{7<t5 


® /)(7CT (g-f) 




-0.538 
-0.230 


-0.461 
0.153 


8 


{7<t5 


(8)€5)(7<T fXi^s) 


-7 

-3 


-5.253 
-2.251 


-6 -4.502 
2 1.501 


9 


{7<t5 


® Cmi'){7o- <8i Cm^) 




2.342 
-0.780 




10 


{7<t5 


® Cmi'){7<t ® Cmp) 




0.224 
-0.074 




11 


{7<t5 


® Cmct)(7it Cmct) 




3.016 
-1.005 




12 


{7<t5 


® ^mct)(7ct Cj/ct) 




0.224 
-0.074 




13 


{7m"^ 


® Cm){7p<t 5<t5) 




0.174 
0.074 


0.149 
-0.049 


14 


{7mi^ 


® 5p){7p<T ® (.vb) 




-0.174 
-0.074 


-0.149 
0.049 


15 


{7mi' 






0.224 
-0.074 




16 


{7mi' 


® (,cTZ,)(lpcT ® Cm) 




-0.224 
0.074 
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)le 8: Renormalization factors for the gauge non-invariant operators VV = (7^ /) (7/^ ® 
in the Landau gauge. Matrix elements are arranged in the same way as in Table. 0. 



VVi VV2 
mixed operator 7 finite 7 finite 



1 


(7m « 


)/)(7m ®?5) 


9 

-3 


37.711/11.392 
-3.708 


26.585/0.266 




z 


(7m 


' ?5j(.7M "X" J) 




n /loa 

U.4ZO 

-0.142 




3 


(7m 


)5m>')(7m ^^mcts) 




-0.174 
-0.074 


-0.149 
0.049 


4 


(7m 


)€mi')(7m ® ^Mi-S) 




-1.813 
-0.777 


-1.554 
0.518 


5 


(7m « 


)C^p)(7m '^?P<t5) 




0.174 
0.074 


0.149 
-0.049 


6 


(7m « 


)Ci.p)(7m iXiCi^ps) 




-1.289 
-0.552 


-1.104 
0.368 


7 


(7<t5 ' 


X>/)(7<t5 


-7 
-3 


-5.253 
-2.251 


-6 -4.502 
2 1.501 


8 


(7<t5 ' 


X>?5)(7^5 (X)/) 




-0.538 
-0.230 


-0.461 
0.153 


9 


(7<t5 ' 


Xl Cmi')(7<t5 ® Ci^ps) 




0.224 
-0.074 




10 


(7<t5 ' 


X> Cm"^)(7<t5 ® Cmi'S) 




3.016 
-1.005 




11 


(7.T5 ' 


Xl Cm't)(7<t5 Cpcts) 




-0.224 
0.074 




12 


(7<t5 < 


X> Cm<t)(7<t5 <X> Cm<^5) 




2.342 
-0.780 




13 


(7mi' 


8>Cm)(7p<t5 8>Cp) 




0.174 
0.074 


0.149 
-0.049 


14 


(7mi' 


^Cp)(7p<t5 ® Cm) 




-0.174 
-0.074 


-0.149 
0.049 


15 


(7mi' 


8iCi^5)(7pct5 fXl^CTs) 




0.224 
-0.074 




16 


(7mi' 


8> C<t5)(7p<t5 ® Ci^s) 




-0.224 
0.074 
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Table 9: Renormalization factors for the gauge non-invariant operators VA = (7;^®/) (7^5® 
^5) and AV = (7^5 (8) /)(7^ ® ^5) in the Landau gauge. Matrix elements are arranged in the 
same way as in Table. ^. 













VAi 




VA2 




AVx 


AV2 




mixed operator 




7 


finite 


7 


finite 


7 


finite 


7 finite 


1 




^m)(75 «>Cm5) 






1.885 
-0.628 








-0.439 
-0.188 


-0.376 
0.125 


2 




Cm)(75 'XlCas) 






0.067 
—0.022 








0.052 
0.022 


0.044 
-0.014 


3 




C<t5)(75 ®Cm) 






0.485 
0.208 




0.416 
-0.138 




0.624 
-0.208 




4 




C<t5)(75 






-3.675 
-1.575 




-3.150 
1.050 




4.853 
-1.617 




5 


{75 5 


Hm)(^'X>Cm5) 






0.565 
-0.188 








-1.466 
-0.628 


-1.257 
0.419 


6 
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OsF 




a, Cir/a (a) b, Dir/a (b) 





Figure 1: One-loop diagrams for bilinear operators. Diagrams (b) and (e) are absent for 
gauge non-invariant operators. 
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b, Di^/a a, Cir/a 




a', C'ir/a (a) 6', D V/" (b) (c) 





(f) 

Figure 2: One-loop diagrams for color two-loop four-quark operators Q"'{js ® ^f)Q°' " 
Q^ilS' ^^F')Q''- Thick horizontal bars at the four-quark vertices signify contraction of spin- 
flavor quantum numbers, while dotted lines represent link factors and flow of color indices. 
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b, Dir/a a, Cir/a 




. (h) 

(f) 

Figure 3: One-loop diagrams for color one-loop four-quark operators Q"'{'ys^^F)Q^-Q^{ls'® 
^F')Q°'- The meaning of the symbols is the same in Fig. |^. 
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